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Abstract 

A numerical integration is made, starting with a Kantowski-Sachs 
(KS) model with an added initial perturbation at time one-tenth the 
age of the universe, and evolving up to the present. According to a 
prediction by Barrow (1989), the ratio between the fluctuations of the 
average matter density and those of the metric tensor varies with the 
inverse squared wavelength of the fluctuation. In this paper we take 
an extravantly large value for this wavelength (20% of the horizon’s 
radius), to see that even then that ratio remains greater than unity. 
So we may be reasonably confident of models assuming density ho¬ 
mogeneity, despite the fact that the scale of observed structure in the 
universe has been increasing lately. 


1 Introduction 

This paper is a sequel to a study on the evolution of small perturbations 
on cosmological models, either of the Friedmann-Lemaitre-Robertson-Walker 


1 


(FLRW) family or closely related to it. The leading idea is to show that 
small perturbations of the spacetime metric are compatible with not-so-small 
perturbations of the average matter density in the universe - cf. Fagundes 
& Kwok (1991). Here we start with a Kantowski-Sachs (KS) model, which 
is similar to the spherical FLRW model, and is easier than the latter to 
deal with in a numerical study: while the space section for the latter is the 
hypersphere S^, for KS it is x R, where 5^ is the ordinary sphere and R is 
the real axis. Then we introduce deviations from the KS metric at an initial 
time t = 0.1, taking the age of the universe as to = 1- The idea is to compare 
this KS model with perturbations which evolve in time with the unperturbed 
KS universe. 

In Fagundes & Kwok (1991), the scale (wavelength) of the perturbation 
was much smaller than the radius of observational horizon, and so, as ex¬ 
pected from Barrow’s (1989) heuristic result, the fluctuations of the metric 
components turned out to be much smaller than the fluctuation of the matter 
density. Here we choose a perturbation with a very large wavelength (1/5 
of the horizon’s radius), with the result that one of the metric components 
still varies less than (about one-half as much as) the density. This result 
reinforces our current conhdence in models based on the assumption of ho¬ 
mogeneity of matter distribution, such as those of the FLRW family, since 
the scale of observed inhomogeneity is believed to be smaller than that of 
the perturbation assumed in this investigation. 

2 Einstein equations 

We choose units such that c = G = to = 1, and assume a metric in the form 

— a^(t,r)dr^ — b‘^{t,r) sin^ r — c^{t,r)dC,‘^ , (1) 

with 0<r<7r,0<(p< 2n, —oo < C < +cxo, which is to be compared with 
the KS metric (Kantowski & Sachs 1966) 

ds^ = dt^ — a%g{t){dr‘^ + sin^ r dip‘^) — c\s{^)dC ■ (2) 

For the metric given by Eq. Einstein equations are (a dot means d/dt, 
a prime d/dr): 

db dc be 1 a'y a'c' b'c' a' 2b' c' b" d' 

ab~^ ac~^ bc~^ a? adb ode d^bc ad a?b a?c a'^b a?c ^ 

(3) 


2 




db' 

ab 


b' 

dd 

d ,d 




- + (-- 

- -) cot r = — 

b 

ac 

c a 

b 

b c 

be 

b'd 

d 




—cot r = 0 

b c 

be 

a%c 

a^c 


(4) 

(5) 


a c ac 

—I-1- 

a c ac 


a'd 


a^c 


a'^c 


= 0 


( 6 ) 


a 0 db 1 a'b' a' 2b' 
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where p = p{t,x) is the average (pressureless) matter density and v = v{t,x) 
is the velocity held with respect to the Hubble how, assumed for simplicity 
to have only the x-component. Terms of the order of are discarded. 

As is well known (see, for example, Weinberg 1972, p. 163), Eqs. | and 
H are constraints on the solutions of the second order quasi-linear Eqs. |- 
^ To solve the latter we dehne new variables f = d, g = b, and h = c, 
in order to get a system of six hrst-order (in t) diherential equations for 
the dependent variables a, b, c, /, g, and h. It is convenient also to dehne 
Jks = b,KS, hxs = cks- 

We want our solution to as close as possible to FLRW’s spherical model 
with a dustlike matter distribution, so we take as our comparison KS solution 
one with this property, namely that in Shikin’s (1967) equations (16) and 
(20), with A = D/3 = d, B = E = 0. (See also Fagundes 1982.) The 
functions axs, cks are obtained in terms of the parameter rj (t), which is the 
inverse of t{ri) = d{ri — sinr^), as in the FLRW models. Here 


aKs{t) = a[l - cosr]{t)], 

CKsit) = 3d{2 - g{t) cot[?7(f)/2]}, 

and we choose d = [tt/S — sin(7r/3)]“^, which corresponds to g{l) = 7r/3, 
^FLRW ^ 4/3^ _ [{A/3)aKs{l)/cKs{l)W = 1.1934557651pS"^ 

where is the usual critical density (present density in Einstein-de Sitter’s 
model). 
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3 Numerical integration 

The system was then numerically integrated, in the interval t = 0.1 to t = 1.0. 
(An attempt was made to start at t = 2 x 10“®, the recombination time; but 
the result was not satisfactory.) The initial conditions are 

a(0.1,r) = aii-s(0.1)[l + q; cos(lOr)] 

&(0.1,r) = aii-s'(0.1)[l + q; cos(lOr)] 
c(0.1,r) = ck 5 ( 0 . 1)[1 +/5cos(10r)] 

/(0.1,r) = /if 5 ( 0 . 1 )[l + Q;cos(10r)] 

5f(0.1,r) = /k 5(0.1)[1 + acos(lOr)] 

/i(0.1,r) = /iif 5 ( 0 . 1 )[l +/5cos(10r)] 

To get the desired effects we put a = (5 = 0.02. Note also the large value 
of the initial perturbation comoving wavelength, L = tt/S. Despite this the 
perturbation remains sub-horizon sized: \{t) ~ < 3t for t > 0.1; so we 

need not worry about gauge ambiguities - see, for example, Kolb & Turner 
(1990). 

The region of integration was divided into 1000 parts for the r variable, 
and 9000 parts for f, which means a spatial interval Sr = 10“^7r and an 
integration step St = 10“^. First and second spatial derivatives were obtained 
by the usual rules, with a shift of 0.5(5r to soften the effect of the coordinate 
singularities at r = 0,7r: for example, for each t = 0.1 -|- nSt, 0 < n < 9000, 
putting r = {k + 0.5)(5r, 0 < A; < 1000, we dehned A[k] = a{t,r), B[k] = 
b(t,r), and so on; and took a'(t,r), b"(t,r) respectively as 

AR[k] = {A[kPlus] — A[kMinus])/2Sr, 
bRR[k] = {B[kPlus] - 2B[k] + B[kMinus])/Sr^, 

where 

f (1, 0) for A; = 0 

{kPlus, kMinus) = < (A; -|- 1, A; — 1) for 0 < A; < 1000 . 

[ (1000, 999) for k = 1000 

After each increment t ^ t + St we make a(t,r) —> a(t,r) -1- f(t,r)St, and 
similarly for b, c, /, g, h. As for p(t,r) and v(t,r) they are given by the 
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constraint equations. Stability of the solution was checked by redoing the 
integration with Sa, 5{3 = ±0.002. 

Calculations followed the simple rules, as given, for example, by Smith 
(1978), and were programmed in the C language on an HP 750 workstation. 

4 Results 

The results at t = 1.0 are shown in Table 1, for uniformly spaced val¬ 
ues of t/tt. The entries are arei = ci{l,r)/axsi.^)) Kei = b{l,r)/aKs{^), 
Orel = c(l,r)/cif 5 (l), n(l,r) in km/sec, and = p(l, r)/p^ 5 (l). Except for 
n(l,r), they were htted to the following short Fourier series, where x = t/tt: 

arei{x) = 1.0008 — 0.1276 cos(107ra;) ± 0.0051 cos(207ra;), (8) 

brei{x) = 0.9932 ± 0.0200 cos(107ra;) ± 0.0069 cos(207ra;), (9) 

Crei{x) = 1.0008 ± 0.0202 cos(107ra;) — 0.0008 cos(207ra;), (10) 

= 1.0236 ± 0.2256 cos(107ra;) - 0.0041 cos(207ra;), (11) 

These hts have the approximate periodicity of the initial perturbation, and 
are symmetrical about p = 7r/2. Not counting the suspicious values near 
r = 0, TT, note the fluctuations of about 13% for arei, 2% for hrei and Creu and 
25% for p^^i. Further details are given in Gonzalez (1994). 

The fluctuation of arei is one-half that of Pj,^i, qualitatively conhrming 
Barrow’s (1989) result, that Sp/p^ ~ {cto/L)^S(f)/(l)Q, where 0 is the Newto¬ 
nian potential (which corresponds to the metric tensor in general relativity), 
and L is the scale of the perturbation. Here we have L = (7r/15) x 3cto, 
or about 20% of the horizon’s radius; this makes Barrow’s estimate Sp/p^ 
~ 2.53 X (50/0Q. In Fagundes & Kwok (1991) we had an initial perturba¬ 
tion with wavelength L = 300h“^ Mpc, and a ratio of the order of 100 was 
obtained between 5p/p^ and 50/00- 

The fact that 5arei is still smaller than 5p^ej, even with our exaggerated 
asumption for the size of L suggests that as long as the scale of observed 
inhomogeneities is smaller than, say, 600/i“^Mpc, models based on the ho¬ 
mogeneity of matter distribuition are reasonably correct as far as the metric 
tensor is concerned. However, this is not an absolute conclusion. In Fagun¬ 
des & Mendonga da Silveira (1995), where initial conditions were contrived 
to illustrate another problem, we got 5p/po ~ 2.5 x 50/0o for L = 120/i“^ 
Mpc. 
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TABLE 1. Results of the integration for the functions dehned in the 
beginning of Section 4. 


r/vr 

^rel 

brel 

^rel 

Prel 

n(l, r) 

0.0005 

0.87197 

0.87197 

1.02292 

1.56328 

1077.29 

0.0504 

0.99773 

0.90637 

1.00227 

1.15528 

-69.32 

0.1004 

1.13339 

0.98062 

0.97988 

0.79305 

-20.79 

0.1503 

0.99464 

1.03031 

1.00157 

0.96130 

82.13 

0.2003 

0.87825 

1.01964 

1.02014 

1.24580 

5.93 

0.2502 

0.99636 

0.98646 

1.00129 

1.02331 

-69.62 

0.3002 

1.13342 

0.98029 

0.97970 

0.79368 

-7.77 

0.3501 

0.99506 

1.00894 

1.00148 

0.99168 

73.20 

0.4001 

0.87831 

1.01973 

1.02015 

1.24548 

1.47 

0.4500 

0.99572 

0.99899 

1.00138 

1.00534 

-70.74 

0.5000 

1.13343 

0.98022 

0.97969 

0.79376 

0.00 

0.5500 

0.99572 

0.99899 

1.00138 

1.00534 

70.74 

0.5999 

0.87831 

1.01973 

1.02015 

1.24548 

-1.47 

0.6498 

0.99506 

1.00894 

1.00148 

0.99168 

-73.20 

0.6998 

1.13342 

0.98029 

0.97970 

0.79368 

7.77 

0.7498 

0.99636 

0.98646 

1.00129 

1.02331 

69.62 

0.7997 

0.87825 

1.01964 

1.02014 

1.24580 

-5.93 

0.8497 

0.99464 

1.03031 

1.00157 

0.96130 

-82.13 

0.8996 

1.13339 

0.98062 

0.97988 

0.79305 

20.79 

0.9496 

0.99773 

0.90637 

1.00227 

1.15528 

69.32 

0.9995 

0.87197 

0.87197 

1.02292 

1.56328 ■ 

-1077.45 
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